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FRICTIONAL FORCES AND COLLISION FREQUENCIES 
BETWEEN MOVING ION AND NEUTRAL GASES 


Peter Stubbe 

Laboratory for Space Sciences 


ABSTRACT 

Two different terms are used in literature for the frictional 
force between moving ion and neutral gases. It is shown how the 
corresponding collision frequencies have to be defined in order to 
keep both terms valid. An expression is derived for the momentum 
transfer collision frequency v in between ions and neutral particles 
as a function of temperature and the relative flow velocity I v. - vj. 
The result indicates that for most practical applications v {n may 
be considered as independent of | v ; ~ y n I • However, for very high 
velocities, e.g. a convection flow from the magnetosphere into the 
ionosphere, this dependence may be significant. Numerical values 
for v iu are presented. In the case of collisions between ions and 
their parent neutral gases, laboratory values are used for the 
resonant charge exchange cross section obtained for high energies, 
but a correction is made for the much smaller thermal energies in 
the temperature range typical for ionospheric conditions. This 
gives rise to an enhancement ofthe collision frequencies, compared 
with the values obtained from the uncorrected charge exchange 
cross sections, the amount of which depends on the particle mass. 
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FRICTIONAL FORCES AND COLLISION FREQUENCIES 
BETWEEN MOVING ION AND NEUTRAL GASES 


1. INTRODUCTION 

The dynamic behaviour of the ionospheric plasma and atmospheric neutral 
gas can be described by the Navier-Stokes equation of hydrodynamics. This 
equation expresses the acceleration of a fluid in terms of a pressure gradient 
force, a viscous force and the resultant of all external forces. One of the most 
important external forces is the frictional force between an ion gas and a 
neutral gas moving with different velocities. The frictional force can be ex- 
pressed in terms of the relative velocity between the ion and neutral gas and 
the number of collisions between ions and neutral particles. 

Two different expressions are used in literature for the frictional force, 
both relating the frictional force to the collision frequency between ions and 
neutral particles. It is necessary, therefore, to define different collision fre- 
quencies in order to keep both expressions valid. One of the goals of the 
present work shall be to give the proper definitions for these collision frequencies. 

Collision frequencies between ions and neutral particles are usually assumed 
to be independent of the relative flow velocity. It shall be investigated in this 
paper to which extent this assumption is correct. 

The frequency of collisions between ions and their parent neutral particles 
is greatly influenced by charge exchange processes. Charge exchange cross 
sections for all processes of atmospheric interest have been measured for 
energies far above the thermal energy range, and they have been extrapolated 
to thermal energies without any correction. In the present paper an attempt 
shall be made to correct the charge exchange cross sections for small energies 
which must lead to an enhancement of the collision frequencies. 


2. EQUATION OF MOTION FOR AN ION GAS 

The Navier-Stokes equation, applied to an ion gas in the earth's atmosphere, 
reads as follows: 


d Fi /l \ 

P. dt - - grad Pj + 7] i U graddivv. + AyJ + k (1) 
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P. = mass density of the ion gas 

p. = partial pressure of the ion gas 

v. = ion velocity 

Tj . = coefficient of viscosity 

k = resultant external force per unit volume. 

The external force k may be written as the sum of a gravity force, an 
electromagnetic force, a Coriolis force, a centrifugal force and a frictional 
force (ion drag force): 

k - P. g + N ; e|E + Vj xBj + 2 / 3 ; |v. xfij + P ; 0 x (rxO) + k ( (2) 

g = acceleration due to gravity 
PL = ion number density 
e = ionic charge 
E = electric field strength 
B = magnetic induction 
R = angular velocity of the earth 

_r = position vector, measured from the earth's center 
k f = external frictional force or ion drag force. 

For a practical application of Equations (1) and (2) to the earth's ionosphere 
the acceleration force, the internal friction force, the Coriolis force, and the 
centrifugal force may be omitted. 


3. FRICTIONAL FORCE AND COLLISION FREQUENCY 

Two different terms for k f are in use in literature. Some authors (e.g. 
Dougherty, 1961; Kendall and Pickering, 1967; Kohl and King, 1967; Priester, 
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Roemer and Volland, 1967) use the expression 


k f 


= -N. v. 


m. 




Pi V; 


hi -sj 


(3) 


while others (e.g. Schluter, 1950; Chandra, 1964) use the expression 


k , - - N . v . /a- fv.-v) (4) 

— f i in 'in ^ — i — n) ' ' 

u. =m.m/m.+m = reduced mass 
m. = ion mass 
m n = neutral particle mass 

v in = number of collisions per second of one ion with the neutral particles 
v n = neutral gas velocity. 

Since v in is a frictitious quantity which merely has the dimension of a 
frequency, it can be defined so that either (3) or (4) are valid. 

Two different definitions for v. n are given in literature. According to 
Banks (1966a), the momentum transfer collision frequency is defined as 



/ hn \ 3 2 r 

\ kT/ 

\2 77 kTj ] 

•'O 


Q(g) 


'2kT g‘ 


g 5 dg 


(5) 


1 'GO 

(1 - cos y)b db = momentum transfer cross section 

o 

g = relative thermal velocity 
X = angle of deflection (see Figures 12 and 13) 
b = collision parameter (see Figure 14) 
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In order to find out how the frictional force term reads when we use v in as 
defined by Equation (5), we refer to the momentum transfer equation provided 
by statistical plasma mechanics. After Burgers (1960), the frictional force in 
the terminology of this paper is given by 


kf = 




477 

T ^in 


N i N n 





Q(g) 


l . /2kT g 2 c 
in/ K 5 


g dg (6) 


or, inserting (5): 


Li. (v. - v ) 

i in r in \ -l n/ 


(7) 


Equation (7) is identical with (4). 

Dalgarno (1961), on the other hand, defines the collision frequency as 


kT 

m i D i 


( 8 ) 


or, when we insert for the diffusion coefficient D. the expression given by 
Chapman and Cowling (1961): 


477 

1 


~~T ' N_ 


^in 

m. 


(w) (skf ( 


Q(g) 


■ n / 


2kT g" 


g 5 dg (8a) 


By comparison of (8a) with (6) we obtain 


k, - -N.z^. m . /v. - v ) 

— f i in i y -l — n ■/ 

Equation (9) is identical with (3). 


( 9 ) 
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The numerical values for v. n , given by Cowling (1945) and Dalgarno (1961 
and 1964), are based on the defining Equations (8) or (8a). Hence, when these 
collision frequencies are used, Equation (3) must be applied for k f . On the 
other hand, the collision frequencies presented by Banks (1966b) and later on in 
this paper (paragraphs 5 and 6), are calculated by means of Equation (5). There- 
fore, these values have to be taken in conjunction with Equation (4). Thereby, 
the apparent discrepancy between the collision frequencies of Dalgarno (1964) 
on the one side and of Banks (1966b) on the other side is removed. 

The frictional force as given by Equation (6) leads, when decomposed into 
the constituents of Equation (4), to an expression for the collision frequency 
which does not depend on the relative flow velocity | v . - vj . Such a dependency, 
however, should be expected, since the velocity distribution functions of the ion 
and neutral gas cannot only be influenced by a change in T, but also by a change 
in ly. - vj . 

We will derive here an expression for k f leading to a collision frequency 
which shows a dependence on |v. -v n | . Our basic assumption is that the ve- 
locity distributions of the ion and neutral gas are Maxwellian distributions, 
displaced by the drift velocities v. and y n , respectively (Smith, 1964): 



Cj, c n = thermal velocities. 

The result of our calculations, which are presented in the appendix, is: 
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By comparison with (4), the momentum transfer collision frequency is 
given by: 



v /2 

^2n kT R J e 


^in/ 2kT R 




e -M in /2l<TR g 2 


g 5 dg (11) 


For v. = y n and T i = T n Equation (11) is identical to (5). We will discuss 
this formula later (paragraph 5). In the following we must distinguish between 
ions in unlike gases and ions in their parent gases. In the first case, the 
momentum transfer cross section is determined by electrostatic interactions, 
while in the second case charge exchange plays a predominant role. 

4. MOMENTUM TRANSFER CROSS SECTION FOR IONS MOVING 
IN THEIR PARENT NEUTRAL GAS 

The momentum transfer section Q was defined as 


Q(g) - 2tt I ( 1 - cos x) b db - 4 tt 
•'0 

X, d and b are defined in Figures 12, 13 and 14. 
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Due to electrostatic interaction forces, which are described later, the in- 
tegrand is greater than zero for all finite values of b and g, but, for given g and 
increasing b, it converges fast enough to zero to guarantee the convergence of 
the above integrals. 

When the ion gas is embedded in its parent neutral gas, the value of Q is 
much less determined by electrostatic interactions than by charge exchange 
processes of the type 


X + + Y ~ X + Y + 


where X and Y are particles of the same neutral gas species. After Holstein 
(1952), the effect of charge exchange processes on the momentum transfer cross 
section can be described as follows: 

When the particle trajectories are linear, a critical collision parameter b c 
with the following features can be defined: For b < b c the charge exchange 
probability P ex is a rapidly changing function of b, oscillating between 0 and 1 
and having an average of 1/2. For b > b c , P ex decreases rapidly with increasing 
b. Therefore, on the average, every second collision with b < b c leads to a 
charge exchange process. Let us consider a particular impact, in the first 
instance without charge exchange (Figure la), then with charge exchange (Fig- 
ure lb). In the first case the ion is deflected by Xi, 

Figures la and lb 


in the second case (since X and Y are of the same species) by \ 2 t'- v - x l (i.e. 
cos x 2 = -cos x t ). Hence, on the average, for b< b c the integrand (1 - cosy) * b 
is equal to b. Therefore, if we neglect the very small contribution of P ox for 
b > b c , we can describe the net result of charge exchange processes in the fol- 
lowing manner: Charge exchange processes between ions and their parent 
neutral particles influence the momentum transfer so as if all collisions with an 
impact parameter b smaller than the critical impact parameter b c occur with a 
deflection angle of 90°. Thus, the momentum transfer cross section can be 
written as 


rrb 2 + 4- 


f 


Q(g) 


2 


co s 2 d b db 


(13) 



The charge exchange cross section S ex , which has been measured for many 
processes of atmospheric interest, is defined as the integral 


S 


e x 



(14) 


where dq is a differential cross section and P ex is the charge exchange probability 
related to dq. Since P e x is a function of b only, S e x can be written as 


S 


e x 


277 



bdb 


(14a) 


Omitting the small contributions of the integrand for b > b c , we get 


S 


e x 




Therefore, Q(g) and S rx are related by 


Q(g) 


2S 


c x 


+ 4 77 



cos 2 $ b db 


(15) 


(16) 


Laboratory measurements (e.g. Stebbings, Smith and Ehrhardt, 1964) as well 
as theoretical studies (Knof, Mason and Vanderslice, 1964) showed that for 
energies above 1 eV, when the trajectories are approximately linear, S ex and the 
relative kinetic energy E -(,u in / 2)g 2 are connected by 

= (C,-C 2 - log 10 E) 2 , (17) 
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where Cj and C 2 are constants. In order to evaluate collision frequencies or 
diffusion constants for ions (e.g. Banks, 1966b and Knof et al., 1964), the authors 
used the unmodified Equation (17). However, at thermal energies characteristic 
of atmospheric conditions, the curvature of the particle trajectories due to 
long range attractive forces should have a marked influence on S ex and Q. We 
will therefore make an attempt to correct the charge exchange cross sections 
for small energies. The corrected values of the charge exchange cross section 
and the critical collision parameter shall be denoted by S e ' x and b c ', respectively, 
to distinguish them from the uncorrected values, given by (15) and (17). 

While the task to relate b c ' to b c for a given relative velocity g can be ex- 
actly solved, an exact relationship between S p ' x and b c ' can be found only when 
all branches of the charge exchange interaction potential are known, but even 
then it is a very difficult problem. We will assume that S e ' x and b c ' are related 
in the same manner as S ex and b c : 



This assumption is in agreement with a proposal by Holstein (1952, p. 835). 
Equation (15a) would be correct to the same extent as (15) if the transformation 
from b' to b were strictly linear. The meanings of b' and b are explained in 
Figure 2. 


Figure 2 


When b c ' is known, Q is given by 


Q(g) 


7Tb ' 2 + 

C 


4tt I cos 2 d b ' db ' , 

Jb c' 


(13a) 


or, because of (15a), by 


Q(g) - 2S ' + 4 t 7 I cos 2 db'db' (16a) 
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We are left now with the problem to calculate cos 2 d as a function of b' and 
g and b c ' as a function of b c and g. 

We assume a long range interaction potential between an ion and a neutral 
particle of the form (e.g. McDaniel, 1964) 


U(r) = 



U(r) = +» 




for r > a 


for r < a 


( 18 ) 


The fourth order term is due to dipole interactions, the sixth order term 
to quadrupole interactions, a is the sum of the gas kinetic radii of the particles, 
A can be expressed in terms of the polarizability v of the neutral particle. 


A = 



(19) 


It would be more realistic to describe the repulsive potential by a term of 
the form C/r 12 , but since no value., for C are available, this would only impede 
the calculations without yielding a higher accuracy. The potential given by 
Equation (18) is illustrated in Figure 3. 

Figure 3 

Using the symbols defined by Figure 4 and applying the conservation laws 

Figure 4 


for energy and angular momentum, we obtain the relations 


u . 

i n 


All n 

+ U(r) - 2 g 

(20a) 

in r 2 j> = M in gb' 

(20b) 
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NJ T3 


w can be expressed as 


w 2 - r 2 + r 2 0 2 


r and <p stand for the time derivatives dr/dt and db/ dt , respectively. Con- 
sidering the identity 


. do 
$ = r d? 


r and do/ dr are given by 


. - 2 > 


1 - 


2U(r) 


( 21 ) 


dep 

dr 


r- 2 


(>d) 


2U(r) 

^in 


- 1/2 


( 22 ) 


By means of (22) we can immediately obtain an expression for ft: 


. x f / b' 2 2U(r)V 1/2 d 


dr 

2 


(23) 


b is the minimum distance between the two particles for given b' and g. If 
r CA (PCA = point of closest approach) is defined as the greatest nullpoint of 
( r ) , b is given by (see Figure 5) 


b 



for r p c a > ° 

for r pcA < a 


(24a) 
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For sufficiently small values of b' , however, there exists no nullpoint of 

Figure 5 


r 2 ( r ) . In this case, b is simply given by 


b - a 


(24b) 


If b is greater than a, i.e. b - r pcA , the integrand in Equation (34) has a 
singularity for r - b. In order to remove this singularity and to make the range 
of integration finite, we change the variables by using the transformation (Mason 
and Schamp, 1958) 


b 

sin 8 - — 


(25) 


We thus obtain for d: 



The critical collision parameter b 0 ' which has the property that for b' < b 0 ' 
r 2 ( r ) has no nullpoint, is of significance for the determination of b c ' , As we 
read off from Figure 5, b ' is given by the condition 


P C A 


( b c') 


(27 a) 


if b c > r P cA (V) andb y 


b 


c 


(27 b) 


if b c < r pcA (b 0 ‘) , Now all formulae necessary for a calculation of Q(g) and 
S e ' x (g) have been provided. Numerical results are presented in the next 
paragraph. 
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5. NUMERICAL VALUES FOR MOMENTUM TRANSFER CROSS 

SECTIONS AND COLLISION FREQUENCIES FOR IONS IN THEIR 
PARENT NEUTRAL GASES 

Experimental values of S ex (g) are available for all processes of at- 
mospheric interest. We will base our calculations on the values for S ex and 
rj given in Table 1. The values for 17 have been taken from Banks (1966b). 


Table 1 

Charge Exchange Cross Sections and Polarizabilities of Neutral 
Particles. Sources: (1) Stebbings et al. (1964); (2) Fite et al. 
(1962); (3) Cramer and Simons (1957); (4) Knof et al. (1964); 

(5) Amme and Utterback (1964). 


Process 

r) (10 24 cm 3 ) 

S ex /2 ( 10 " 8 Cm ) 

Source for S ex 

0 + + 0 - 0 + o + 

0.89 

5.95 - 0.63 log 1Q E(eV) 

(1) 

H + + H - H + H + 

0.67 

7.60 - 1.06 log j Q E(eV) 

(2) 

He + He - He + He + 

0.21 

5.25 - 0.74 log j 0 E(eV) 

(3) 

N + + N - N + N + 

1.13 

5.53 - 0.46 log j 0 E(eV) 

(4) 

o 2 + + o 2 — o 2 + o 2 + 

1.60 

5.37 - 0.54 log j Q E(eV) 

(5) 

n 2 + + n 2 - n 2 + n 2 + 

1 

1.76 

7.36 - 0.68 log j 0 E(eV) 

(5) 


Using the formulae derived in paragraph 4, Q(g) and S e ' x (g) have been 
calculated numerically with the help of an electronic computer. Since B (de- 
fined by Equation (18)) is practically unknown, a set of values for B, ranging 
from 0.5 to 2.0 10 " 8 cm has been used. The result was that the influence of 
the sixth order potential term is completely negligible, since in no case was 
the deviation in Q(g) caused by this term more than 0.5%. Therefore, all 
further calculations have been carried out with B = 0, i.e. for a pure fourth 
order potential law. 

Results for Q(g) and S e ' x (g), compared with S ex (g), for the processes 
0 + + O - O + Q + and H + + H - H + H 4 are shown in Figures 6 and 7. 
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Figure 6 
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Figure 7 

The numerical values for Q(g) have been used to claculate the momentum 
transfer collision frequency v in after Equation (11) for the processes listed in 
Table 1. Figure 8 shows v in for the most important charge exchange collision 
process in the ionosphere, namely 0 + + O -* O + 0 + , as a function of the re- 
duced temperature T R for | y £ — y n S =0. In order to indicate the effect of 
electrostatic interactions, v in has also been calculated for the uncorrected 
momentum transfer cross section Q = 2 S e x , which would be the case if no 
electrostatic interactions would be present. 

Figure 8 


Since the difference between Q and 2S ex increases with decreasing velocity, 
the influence of electrostatic interactions is conspicuous only for low tempera- 
tures (below 1000 °K), and is more noticeable for heavier particles because the 
mean thermal velocity is proportional to m~ 1/2 . For the process 0 2 + + O 2 ^ 0 2 + 0 2 + , 
for instance, v in , calculated with the exact Q(g) , is greater than v in for Q = 2S CX 
by 39% for 400 °K and 17% for 1000 °K. On the other hand, the corresponding num- 
bers for the process H + +H- > H + H + are 2.3% and 0.4%. 


In the temperature range 500°K < T R < 3000°K the collision frequency v i n for 
| v. - v n | = 0, as a function of T R , can be approximated by the following expres- 
sions, when the neutral particle number density is measured in the unit cm' 3 : 


sec 


v( 0 + , O) = 1.86 • 10~ 9 (T R /1000) 0 37 n(O) 

v(H + , H) = 12.03 • 10' 9 (T r / 1000)°- 38 n(H) sec -1 

v (He + , He) = 2.92 • 10' 9 (T R /1000)°- 37 n(He) sec' 1 

v (N + , N) = 1.75 • 10~ 9 (T R /1000) 0-34 n(N) sec' 1 

^(0 2 + , 0 2 ) = 1.17 • 10' 9 (T R /1000)° 28 n(0 2 ) sec' 1 

z^(N 2 + , N 2 ) = 2.11 • 10~ 9 (T r / 1000)°- 38 n(N 2 ) sec' 1 

These values can be used to calculate the ambipolar diffusion constant D a for 
an electrically neutral electron-ion gas in the ionosphere, which is given by 


D. 


k (Tj + T e ) 

v . a. 

in / in 


sin 2 I 
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In the F2 region, where 0 + and O are the major ionic and neutral constituents, 
D a can be written as 


D 


= 0.55 • 10 16 


1000 


0 . 37 


(T, + T.) 


sin 


n(O) 


2 - 1 
cm sec 


In the exosphere, where H + and H are predominant, D a is 


D 

a 


1.38 • 10 16 



(l i + T e) Sin 2 I 

n(H) 


cm sec 


v in as a function of | v. - v n | for the processes 0 + + O - O + 0 + and H + + H 
- H + H + is shown in Figures 9 and 10. When we consider diffusional flows, 
electromagnetic drifts or neutral gas winds in the upper atmosphere and 
ionosphere, | v L - v | should not exceed about 300 m/sec. 

Figure 9 


Figure 10 

Therefore, when applied to these mechanisms, v i n may be considered as 
independent of |y. - v n |. However, for particle flows from the magnetosphere 
into the ionosphere, as proposed by Axford and Hines (1961), this dependence 
may be of significance, since the plasma convection velocity is expected to be 
of the order of kilometers per second. 


6. NUMERICAL VALUES OF COLLISION FREQUENCIES FOR IONS 
IN UNLIKE NEUTRAL GASES 

For reasons of completeness, collision frequencies for ions in unlike neutral 
gases are also presented, although the basic studies in this field have been car- 
ried out more than six decades ago by Langevin (1905) and later on by Hasse 
(1926). Langevin and Hasse" give an expression for the ion mobility K which, 
using the terminology of this paper, can be written as 


K 


A(M 

2N n n ^ 


(28) 
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where A(A) , A being defined as 


A. 



(29) 


is given by Hasse (1926, Table III). A is not identical with the quantity defined 
by Equation (18). Using Equations (1), (2), (4) and taking v n = 0, K andv^ are 
found to be related by 


2 * F(A) 



(30) 


F(A) - /rr/A(A) is shown in Figure 11. 

Figure 11 


Since a, the sum of the gas kinetic particle radii, is not known well enough, 
the usual method is to take the limiting value F(0), i.e. for vanishing a or in- 
finite 77, instead of the correct value F(A). This leads to the approximation 


v . 

in 



(31) 


We must realize, however, that the values obtained from Equation (31) can 
easily be wrong by about 20% for temperatures below 2000°K and by more than 
this for higher temperatures. Especially the temperature dependence of v in is 
not correctly described by the approximation (31). For A < 0.6 the collision fre- 
quencies decrease with increasing temperature, while for A > 0.6 the opposite is 
true. For A > 1.2 the collision frequencies approximately follow a T 1 2 -law. 
Assuming T = 500°K and 77 = 10' 24 cm 3 , a value of A A 1.2 is adopted for a^. 4 
10' 8 cm. iA n after Equation (31) is based on the defining Equation (5). The 
defining Equation (8a) would yield 


v . 

1 n 


= 6.94 



N 


(31a) 
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Although momentum transfer collision frequencies can simply be evaluated 
from Equation (31), we will give some numerical values for the most important 
collision processes in the ionosphere: 


^( 0 + . 0 2 ) = 1.00 
v (O + , N 2 ) = 1.08 
v(0 + , H) = 2.19 
v(H + , O) = 2.52 
v (NO + , 0 2 ) = 0.83 
i'(NO + , N 2 ) = 0.90 
v(NO + , O) = 0.76 
^(0 2 + , N 2 ) = 0.89 
v (O 2 + , O ) = 0.75 


10' 9 n(0 2 ) 

sec -1 

10' 9 n(N 2 ) 

sec -1 

10 ' 9 n(H) 

sec -1 

10' 9 n(O) 

sec -1 

10' 9 n(0 2 ) 

sec -1 

10" 9 n(N 2 ) 

sec -1 

10~ 9 n(O) 

sec -1 

1(T 9 n(N 2 ) 

sec -1 

10" 9 n(O) 

sec -1 


As mentioned before, the collision frequencies defined by (5) and (8a) are 
fictitious quantities, that is, they are quantities having the dimension of a fre- 
quency and were introduced merely for reasons of convenience. When two 
particles approach each other so closely that they are separated only by short 
range repulsive forces or, in other words, when two particles have a direct 
contact, we will call this a "real" collision. According to this definition, the 
number of real collisions per second, v R , is given by (Stubbe, 1966) 



where u Q and v 0 are abbreviations standing for the expressions 


U 


0 
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-1 L2 

cr 2 U 


in 


1/2 


In terms of A., defined by (29), can be written as 


= HZ 


. 1/2 


r iA 

N 

2 

1 2 A. ,,\2 

e u du + e 1 

/ " 


i 

o 


(32a) 


For temperatures below about 2000°K, (32) and (32a) can be approximated by 


V 


R 



(32b) 


By comparison of (31), (31a) and (32b) we see that v R lies between v. m after 
(31a) and u in after (31). 

8. SUMMARY 

a) The relationship between the frictional force and the collision frequency 
between ion and neutral gases is studied. Proper definitions are given for the 
collision frequency to fit the expressions used for the frictional force. 

b) An expression is derived for the momentum transfer collision frequency 
v. as a function of the relative flow velocity |v. - v I . It is shown that for most 
practical applications v in may be considered as independent of |v. -y n l , but 
that for velocities in the km/sec range, as they may be expected in convection 
flows from the magnetosphere, this dependence is significant. 

c) Numerical values of v. for the most important collision processes in 
the ionosphere are presented. In the case of ions in their parent neutral gases, 
resonant charge exchange is of great influence on the collision frequencies. 
Laboratory values are used for the charge exchange cross sections, although 
obtained for higher energies, but an approximate correction is made for thermal 
energies in order to take into account the effect of curved particle trajectories 
on the charge exchange cross section. 
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APPENDIX 

DERIVATION OF EQUATION (10) 

The frictional force between two gases is the result of the momentum trans- 
fer between the colliding particles. We introduce the following quantities: 

c it c = ion and neutral particle velocity before the collision 

c , c r ' = ion and neutral particle velocity after the collision 

g = c . - relative velocity before the collision 

g' = c relative velocity after the collision 

I = m.c. ion momentum before the collision 

— 1 l — 1 

I = m. c.' ion momentum after the collision 

— l 1—1 

AI = I ' -I. momentum transfer 

— in — l — l 

Using the conservation law for momentum, it can easily be derived that 
AI in is given by 

A Iin = ^in U' ( 37 ) 

The conservation law of energy which holds for elastic collisions gives 

Igl = Ig'l = g (38) 

In order to determine g' -g we will distinguish between direct collisions 
(direct contact of the colliding particles) and indirect collisions (no direct con- 
tact of the particles, but curved orbits due to long range electrostatic forces). 

a) Direct Collisions: Under the assumption that no angular momentum is trans- 
ferred during the collision, the orbit is symmetrical to the radial unit 


21 



Figure 12 


vector e r . According to Figure 12, we therefore obtain: 

£ ' _ J S = I £* _ gl ’ e r 

Is' ~ k \ ~ 2g cos d 

b) Indirect Collisions: According to Figure 13: 

£' ~ £ = “ Is' "sl ' e r 

lg' -gl = “ 2g cos d 

Figure 13 

In both cases we get the same result, namely 

s' " S = 2g cos d e r (39) 

or 

A I in = 2/x. n gcosde r (40) 

We now introduce a Cartesian coordinate system (Figure 14) and we es- 
pecially assume that the relative macroscopic velocity v ; - v has the direction 
of the negative z-axis. In this case we are interested only in the z-component 

Figure 14 

of AF n since, on the average, all other components cancel out because of the 
complete symmetry about the z-axis. 

( Al in) z = 2/^in g cosd e rz (41) 
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t 


It can easily be shown that 


e fz - cos $ cos 6 + cos a sin d sin 9 


Hence, the z-component of the momentum transferred from an ion to a neutral 
particle for any particular collision is given by 

K„), = 2K, n g ( cos 2 d cos 9 +cos ttsin i?cos d sin 9 ) (42) 

(AI in ) z and the frictional force k f are connected by 


dk , - ( AI . ) N . dv . e 

— f \ ln/zi l — z 


(43) 


dv. is the number of collisions of one ion per second impinging on the differential 
cross section dq = b • db • da (see Figure 14) and occurring in the velocity range 
g • • • g +dg and the angular range 9 • • • 6 + d<9. According to this definition, dv . is 
given as the product of the relative velocity g, the cross section dq, and the num- 
ber of neutral particles dN n per unit volume having a velocity between g and g + dg 
and being in the angular range between 9 and 9 + 69 , related to the particular ion 
under consideration. 

In order to determine dN n , we assume that the ion gas and the neutral gas 
have Maxwellian velocity distributions, displaced by the velocities y t andv n , 
respectively. Furthermore we assume, in the first instance, that the ion tem- 
perature and the neutral gas temperature are equal. 


dN. 

1 

"N~ 

1 



(44a) 


dN 

n 

"n~ 


n 



e 


~ m r/ 2kT ( ■£, 


) 2 


d 3 


c 

— n 


(44b) 
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K 


We multiply the distribution functions with each other and introduce a new 
velocity 


U 


m . c . + m c 

i — i n — n 

m . + m 

i n 


After some simple manipulations we get: 


dN . dV = N. 


/ m i m n V 2 -l/2kT[(m i + mn )(U-m i v i + m i ^/m i +m n ) 2 +M in (y i -v o -g) 2 ] , 3 

1 n \4^ 2 k 2 TV 


£ 


The integration over U can simply be carried out yielding 


dN. dN 

i n 


N. N 


2rrkT 


3/2 


- n /2kT(v-v n - E ) 2 d 3 


£ 


(45) 


When we drop the restriction T T. = T n we obtain the same result provided 
that we replace T with the reduced temperature T R defined by 



(46) 


Using Figure 14 and defining 


X 


Mi 


kT r 


v . - v 
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we obtain 


e "^in/ 2kT R(^rin'6) 2 = e "^ in/ 2kT R J 2 _ e "M in /2kT R g 2 _ ^ in /kT R g | V .~v „ 1 cosS 


^in/ kT Rgl^i'Inl cos0 


1 + X g cos 0 + yi X 2 g 2 cos 2 Q + yr X 3 g 3 cos 3 0 + 


d 3 g - g 2 sin 9 66 d<£dg 


Inserting these expressions into (45), carrying out the integration over 
<p (0 ^ 4 > — 2 tt) and dividing by Pb , we get for dN n : 


\ 3/2 II2 : 

■ 1 -^ in / 2kT R v i** v n 2 “M in /2kT g 

dN_ = N_ • 2tt ) -e in/ R - 1 - n • g 2 • e ir,/ 


n n 


^sin 8 + X g sin 6 cos 6 + 27 X 2 g 2 sin 8 cos 2 6 + • • dg d8 (47) 


du is given by 


dv . - p * b * dN * db * da 

1 » n 


(48) 


Hence, after integrating over a(0 - 2 tt): 


dk f = e z 4v(± in N. g 2 b cos 2 d cos 6 dN n db 


(49) 
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d is a function of the collision parameter b and the relative velocity g. We 
introduce a velocity dependent collision cross section Q(g) defined by 


Q(g) 



cos 2 d b db 


( 50 ) 


where Q(g) is identical to the momentum transfer cross section 


Q(g) 



( 1 - cos x) b db 


(50a) 


commonly used in literature since, according to Figures 12 and 13, the deflec- 
tion angle x is related to d by means of 


2 cos 2 d - 1 - cos x 


We integrate over 6(0^6 ^-n) and over b(0 ^b <®) in Equation (49), insert 
Q(g) given by (50), consider the relationship 


e 

— z 


V . 
— i 


“ V 
— n 




“ V 
— n 


and integrate over g(0 ^-g < oo) to get the final result 
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FIGURE CAPTIONS 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


la, b. Left side (la): Collision without charge exchange 
Right side (lb): Collision with charge exchange 

2. Illustration of some quantities used in the text 

3. Interaction potential U as a function of the particle distance r 

4. Illustration of some quantities used in the text 

5. r 2 (r = radial velocity) as a function of the particle distance r 

for three different collision parameters b' 

6. Q, 2S e ' x and 2S ex as a function of the relative velocity g for the 
process 0 + + O - O + 0 + 

7. Q, 2S ex and 2S ex as a function of the relative velocity g for the 
process H + + H* - H + H + 

8. v(0 + , 0)/n(0) as a function of the reduced temperature T R for 
lYj - v n | =0. Upper curve for the momentum transfer cross 
section Q after Equation (16a), lower curve for the uncorrected 
momentum transfer cross section 2S 

ex 

9. v(0 + , 0)/n(0) as a function of \v i -y n I for three different 
temperatures 

10. y(H + , H)/n(H) as a function of | y i - Y n I for four different 
temperatures 

11. F as a function of A. 

12. Relation between g' - g and# for direct collisions 

13. Relation between g' - g and d for indirect collisions 

14. Illustration of some quantities used in the text 
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ure 5- (r = radial velocity) as a Function of the Particle Distance r for Three Different Collision Parameters 
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Figure 11. F as a Function of X. 
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Figure 12. Relation Between ' -g_ andi?for Direct Collisions 
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Figure 14. Illustration of Some Quantities Used in the Text 
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